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Abstract 

Based on quantum field theory at finite temperature we carried out new results for 
two-particle Bose-Einstein correlation (BEC) function C2{Q) in case of light hadrons. 
The important parameters of BEC function related to the size of the emitting source, 
mean multiplicity, stochastic forces range with the particle energy and mass depen- 
dence, and the temperature of the source are obtained for the first time. Not only 
the correlation between identical hadrons are explored but even the off-correlation 
between non-identical particles are proposed. The correlations of two bosons in 4- 
momentum space presented in this paper offer useful and instructive complimentary 
viewpoints to theoretical and experimental works in multiparticle femtoscopy and in- 
terferometry measurements at hadron colliders. This paper is the first one to the 
next opening series of works concerning the searching of BEC with experimental data 
where the parameters above mentioned will be retrieved. 



1 Introduction 



For the aim to explore the correlations of Bose-Einstein type (BEC) one needs to use 
the properties of a particle detector, e.g., its tracking system to study the hadron 
processes at some energy region. Such a study will be done soon in the next papers. 

This paper describes an attempt to address the problems of BEC within the the- 
oretical aspects prior the real data will be analyzed. 

Over the past few decades, a considerable number of studies have been done on the 
phenomena of multi-particle correlations observed in high energy particle collisions 
(see the review in [1]). It is well understood that the studies of correlations between 
produced particles, the effects of coherence and chaoticity, an estimation of particle 
emitting source size and the temperature play an important role in this branch of 
high energy physics. 

By studying the Bose-Einstein correlations of identical particles (e.g., like-sign 
charge particles of the same sort) or even off-correlations with respect to different- 
charge bosons, it is possible to predict and even experimentally determine the time 
and spatial region over which particles do not have the interactions. Such a surface 
is called as decoupling one. In fact, for an evolving system such as, e.g., pp collisions, 
it is not really a surface, since at each time there is a spread out surface due to 
fluctuations in the final interactions, and the shape of this surface evolve even in 
time. The particle source is not approximately constant because of energy- momentum 
conservation constraint. 

More than half a century ago Hanbury-Brown and Twiss [2] used BEC between 
photons to measure the size of distant stars. In the papers in [3] and [4] , the master 
equations for evolution of thermodynamic system created at the final state of the 
(very) high multiplicity process were established. The equations have the form of 
the field operator evolution equation (Langevin-like [5]) that allows one to gain the 
basic features of the emitting source space-time structure. In particular, it has been 
conjectured and further confirmed that the BEC is strongly affected by non-classical 
off-shell effect. 

The shapes of BEC function were experimentally established in the LEP experi- 
ments ALEPH [6], DELPHI [7] and OPAL [8], and ZEUS Collaboration at HERA [9], 
which also indicated a dependence of the measured so-called correlation radius on the 
hadron (tt, K) mass. The results for tt^tt^ and tt^tt^ correlations with pp collisions 
at = 1.8 TeV were published by E735 Collaboration in [10]. 

One of the aims of this paper is to carry out the extended model of BEC in the 
framework of quantum field theory at finite temperature {QFTfj) approach which to 
be applied later to real experimental data on two-particle BEC. It is known that the 
effective temperature of the vacuum or the ground state or even the thcrmalizcd state 
of particles distorted by external forces is occurring in models quantized in external 
fields. One of the main parameters of the model considered here is the temperature 
of the particle source under the random source operator influence. 

Among the results obtained in this paper we mention a theoretical estimate ac- 
cessible to experimental measurements of two-particle BEC and proof that quantum- 
statistical evolution of particle- antiparticle correlations are not an artifact of the stan- 
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dard formalism but a quite general properties of particle physics. The effect (called 
as surprized one) for non- identical particles correlations was predicted already in [11]. 



2 Two-particle BEC 

A pair of bosons with the mass m produced incoherently (in ideal nondisturbed, 
noninteracting cases) from an extended source will have an enhanced probability 
C2{pi,P2) — Ni2{pi,p2)/[Ni{pi) ■ N2{p2)] to be measured (in terms of differential 
cross section a), where 

N12{PUP2) - - 1 
(7 dill M"2 

to be found close in 4- momentum space 3^4 when detected simultaneously, as compared 
to if they are detected separately with 

The following relation can be used to retrieve the BEC function C2 (Q) : 

^"(^) = W(i)' ^'^• = +'-'0' 

where N"^^ (Q) in general case refer to the numbers N^^{Q) for like-sign charge 
particles (eg., tt^tt^, K^K^, ...); N^^{Q) — for different charge bosons (eg., 
7r=t7r=F, K'^K'^, ...) or even for neutral charge particles iV"''((3) (eg., 7r°7r°, K^K^, . . .) 
with 



Q = sJ-iPi - P2)m • (Pi - P2Y = v/M2-4m2. (4) 
In formula ([3]) and (jl]) N^'^^ is the number of pairs without BEC and p^^ii = 1, 2) 



are four-momenta of produced particles, M ~ \J {pi +^2)^ is the invariant mass of 

the pair of bosons. For reference sample, N'^'^f [Q), the like-sign pairs from different 
events can be used. It is commonly assumed that the maximum of two-particle BEC 
function C2{Q) is 2 for pi — p2 if no any distortion and final state interactions are 
taking into account. 

In general, the shape of BEC C2{Q) function is model dependent. The most simple 
form of Goldhaber-like parameterization for C2{Q) [12] has been used for data fitting: 

C2{Q) = Co • (1 + Ae-Q'^') • (1 + £Q), (5) 

where Co is the normalization factor, A is so-called the chaoticity strength factor, 
meaning A = 1 for fully chaotic and A = for fully coherent sources; the parameter 
R is interpreted as a radius of the particle source, often called as the "correlation 
radius" , and assumed to be spherical in this parameterization. The linear term in 
dS]) is often supposed to be account for long-range correlations outside the region of 
BEC. However, the origin of these long-range correlations as well as the value of e 
are unknown yet. Note that distribution of, e.g., pions and kaons can be far from 
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isotropic, usually concentrated in narrow jets, and further complicated by the fact 
that the light particles with masses less than 1 GeV often come from decays of long- 
lived heavier resonances and also are under the random chaotic interactions caused 
by other fields in the thermal bath. In the parameterization ^ all of these problems 
are embedded in the random chaoticity parameter A. 

We obtained the C2{Q) function within QFTp approach [3] in the form: 

C2(Q)=eW- [i + ^Y^^ V^+^Y^^f)(Q)] -i^W-Ax), (6) 

where £,{N) depends on the multiplicity N as 

^ '-^^ in 

The consequence of the Bogolyubov's principle of weakening of correlations at large 
distances [13] is given by the function Aa;) of weakening of correlations at large 
spread of relative position Ax 

normalized as F{Q,Ax = oo) ~ 1. Here, f{Q,Ax) is the two-particle distribution 
function with Ax, while f{pi) are one-particle probability functions with i — 1,2; rf 
is a measure of weakening of correlations with Ax: rf as Ax —^ oo. 

The important parameter a in ([H]) summarizes our knowledge of other than space- 
time characteristics of the particle emitting source. 

The fl{Q) in ^ has the following structure in momentum space 

n{Q)^n{Q)-j{n), (9) 

where 

n{Q) = exp(-Ap5R) = exp - pa)^ (pi - P2T] (10) 

is the smearing smooth dimensionless generalized function, is the (nonlocal) 
structure tensor of the space-time size (BEG formation domain), and it defines the 
spherically-like domain of emitted (produced) particles. 

The function 7(71) in Q reflects the quantum features of EEC pattern and is 
defined as 

^ n{uj) = n{uj,l3) ^ -, uj ^ — - — , 11) 

n{ijj) n[Lu') e^'^ — 1 2 

where n{uj,P) is the mean value of quantum numbers for BE statistics particles with 
the energy lo and the chemical potential in the thermal bath with statistical equi- 
librium at the temperature T = The following condition f n f (lu , f3) = N is 
evident, where the discrete index / reflects the one-particle state /. 

Note that it is commonly assumed for a long time that there are no correlation 
effects among nonidentical particles (e.g., among different charged particles). This 
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assumption is often used in normalizing the experimental data on C2 with respect 
to C2"'. In the absence of interference or correlation effects between, e.g., tt"*" and 7r~ 
mesons it is supposed that C^^ = 1. 

In terms of time-like Rq, longitudinal Rl and transverse Rt components of the 
space-time size R^ the distribution A^jj looks like = +, — , 0) 

- ^% = {^p'fRl + (Ap^)^i?i + (Ap^)2i?2 . (12) 

Seeking for simplicity one has {Rl — Rt = R) 

A^fl = (P? - P^fRl + {Pi - P2rR' (13) 
for like-sign charge bosons, while 

^% = {pI + plfRl + ipi + P2rR' (14) 

for different charge particles. 

Obviously, the BEC effect with Q,^^ = exp(— A^^) is smaller than that defined 
by il" = exp(— Apjj). The distribution fi*-' gives rise to an off-correlation pat- 
tern between different charge particles. The evidence of C2 correlation represents 
a quantum-statistical correlation between a particle and an antiparticle. Since we did 
not follow special assumptions on the quantum operator level for C2 from the initial 
stage, it may correspond to a physically real and observable effect. This pattern may 
lead to a new squeezing state of correlation region. We obtain that within the QFTp 
the BEC is more generally sensitive to particle- antiparticle correlations than it would 
be expected from the two-particle (symmetrized) wave function which never leads to 
such the correlations. 

3 Green's function 

In this paper, we would like to focus on the role of the particle mass, which influences 

the correlations between particles. To explore this problem, one wist derive the 
memory history of evolution of particles produced in high energy collisions using the 
general properties of QFT at finite temperature. 

We consider the thermal scalar complex fields $(.t) that correspond to tt^ mesons 
with the standard definition of the Fourier transformed propagator F[G{p)\ 

F[G{p)] = G{x -y) = Tr {r[$(a;)$(t/)]p^} , (15) 

with = e-^"/Tre-^" being the density matrix of a local system in equilibrium 
at temperature T = (3~^ under the Hamiltonian H. 

We consider the interaction of i>(x) with the external scalar field given by the 
potential U. In contrast to an electromagnetic field, this potential is a scalar one, but 
it is not a component of the four- vector. The Lagrangian density can be written 

L{x) = d^,^*{x)d''^x) - (m^ + U)^*{x)^x) 
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and the equation of motion is 



{S/^ +m^)^{x) = -J{x), (16) 

where J{x) — U^{x) is the source density operator. A simple model like this allows 
one to investigate the origin of the unstable state of the thermalized equilibrium 
in a nonhomogeneous external field under the influence of source density operator 
J{x). For example, the source can be considered as (5-like generalized function J{x) — 
ji p{x,e)(^{x) in which p{x,e) is a (5-like succession giving the (5-function as e — > 
(where fl is some massive parameter). This model is useful because the 5-like potential 
U{x) provides the model conditions for restricting the particle emission domain (or 
the deconfinement region). We suggest the following form: 

j(x) = -s(ia^) <p{x) + Jr{x), 

where the source J{x) decomposes into a regular systematic motion part T,{idf^) $(x) 
and the random source Jr{x). Thus, the equation of motion (jl6p becomes 

[V^ + - I](j9^)]$(a:) = -Jr{x), 

and the propagator satisfies the following equation: 

[-pl + v?~t{p^)]G{p^)^l. (17) 

The random noise is introduced with a random operator ri{x) = — Yiiid^), for 
that the equation of motion looks like: 

{V^ + m2[i + j^{x)]}<^{x) = ~Jr{x). (18) 

We assume that rj{x) varies stochastically with the certain correlation function 
(CF), e.g., the Gaussian CF 

(77(a;)r/(j/)> = Cexp(-zVc/i)> z^x~y, 

where C is the strength of the noise described by the distribution function exp(— z^/i^^) 
with Lch being the noise characteristic scale. Both C and define the influence of 
the (Gaussian) noise on the correlations between particles that " feel" an action of an 
environment. The solution of Eq. (fT5|) is 

Hx) = - J dyG{x,y)jR{y), (19) 

where the Green's function obeys the Eq. 

{V^ + m^[l + rj{x)]}G{x, y) = 5{x - y). 

The final aim might having been to find the solution of Eq. p^ . and then average 
it over random operator ri{x). Note that the operator M(x) = V'^ -I- m?[\ -f ri{x)] in 
the causal Green's function 

G{x, y) = . \ . S{x - y) 
M [x) + 10 
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is not definitely positive. However, we shall formulate another approach, where the 
random force influence is introduced on the particle operator level. 

We introduce the general non-Fock representation in the form of the operator 
generalized functions 

b{x) ^ a{x) +r{x), (20) 
b+{x) ^ a+{x) +r+{x), (21) 
where the operators a{x) and a+(a;) obey the canonical commutation relations (CCR): 

[a{x),aix')] = [a+{x),a+{x')] = 0, 

[aix),a+ix')] = S{x-x'). 

The operator-generaHzed functions r{x) and r^{x) in (PO)) and (pij) . respectively, 
include random features describing the action of the external forces. 

Both 6+ and b obviously define the CCR representation. For each function / from 
the space S{dioo) of smooth decreasing functions, one can establish new operators 
6(/) and b+{f) 

Hf) = / f{x)b{x)dx = aif) + [ f{x)r{x)dx, 



b+{f)= / mb+{x)dx = a+if)+ / fix)r+{x)d: 



The transition from the operators a{x) and a'^{x) to b{x) and b'^{x), obeying those 
commutation relations as a{x) and a~^{x), leads to linear canonical representations. 



4 Evolution equation 

Referring to [3] for details, let us recapitulate here the main points of our approach in 
the quantum case: the collision process produces a number of particles, out of which 
we select only one (we assume for simplicity that we are dealing only with identical 
bosons) and describe it by stochastic operators b{p,t) and b~^{p,t), carrying the fea- 
tures of annihilation and creation operators, respectively. The rest of the particles are 
then assumed to form a kind of heat bath, which remains in an equilibrium charac- 
terized by a temperature T (one of our parameters) . We also allow for some external 
(relative to the above heat bath) influence on our system. The time evolution of such 
a system is then assumed to be given by a Langevin-type equation [3] for stochastic 
operator b(p, t) 

idtb{p, t) = A{p, t) + F(p, t) + P (22) 

(and a similar conjugate equation for b'^{p,t)). We assume an asymptotic free undis- 
torted operator a{p,t), and that the deviation from the asymptotic free state is pro- 
vided by the random operator r(p, t): a{p, t) b{p, t) = a{p, t) + r{p, t). This means, 
e.g., that the particle density number (a physical number) {n{p,t))^f^ = (^(p)) + 0(e), 
where (n(p, t))^^^ means the expectation value of a physical state, while {n{p)) denotes 
that of an asymptotic state. If we ignore the deviation from the asymptotic state in 
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equilibrium, we obtain an ideal fluid. One otherwise has to consider the dissipation 
term; this is why we use the Langevin scheme to derive the evolution equation, but 
only on the quantum level. We derive the evolution equation in an integral form that 
reveals the effects of thermalization. 

Equation ([22|) is supposed to model all aspects of the hadronization processes (or 
even deconfinement). The combination A{p, t) + F{p, t) in the r.h.s of (|^^ represents 
the so-called Langevin force and is therefore responsible for the internal dynamics of 
particle emission, as the memory term A causes dissipation and is related to stochastic 
dissipative forces [3] 

/ + 00 
dTKip,t-T)b{p,T) 

with K{p, t) being the kernel operator describing the virtual transitions from one 
(particle) mode to another. At any dependence of the field operator K on the time, 
the function A{p, t) is defined by the behavior of the system at the precedent moments. 
The operator F{p, t) in ([22|) is responsible for the action of a heat bath of absolute 
temperature T on a particle in the heat bath, and under the appropriate circumstances 
is given by 

/ + 00 7 
^ ^V^(p^)c(p,,)e— *. 

The heat bath is represented by an ensemble of coupled oscillators, each described 
by the operator c(p^) such that [c(p^), c"'"(pj^)] ~ ^"^{Pfj. — p'^), and is characterized 
by the noise spectral function ipiPfj.) [3]- Here, the only statistical assumption is that 
the heat bath is canonically distributed. The oscillators are coupled to a particle, 
which is in turn acted upon by an outside force. Finally, the constant term P in (|22p 
(representing an external source term in the Langevin equation) denotes a possible 
influence of some external force. This force would result, e.g., in a strong ordering of 
phases leading therefore to a coherence effect. 

The solution of equation ([221) is given in S{di4) by 

^(p^) = — hr^ [^(^^A^) + pl'^p, e)], (23) 

uj - K{pf,) 

where uj in p{uj, e) was replaced by new scale ujp ^ uj / P . It should be stressed that the 
term containing p{ujp, e) as e ^ yields the general solution to Eq. (|^^ . Notice that 
the distribution p{ujp,e) indicates the continuous character of the spectrum, while 
the arbitrary small quantity e can be defined by the special physical conditions or the 
physical spectra. On the other hand, this p{wp, e) can be understood as temperature- 
dependent succession p{uj, e) = J dxexp{iuj — e)x — > (5(w), in which e —> l3~^. Such a 
succession yields the restriction on the /3-dependent second term in the solution ([^5]) , 
where at small enough T there is a narrow peak at w = 0. 

From the scattering matrix point of view, the solution ((25)) has the following 
physical meaning: at a sufficiently outgoing past and future, the fields described by 
the operators d{pfj) are free and the initial and the final states of the dynamic system 
are thus characterized by constant amplitudes. Both states, </?(— oo) and Lp{+oo), are 
related to one another by an operator S{f) that transforms state Lp{—oo) to state 
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ip{+oo) while depending on the behaviour of f(j)fj): 

ip(+oo) — S{f)ip{—oo). 

In accordance with this definition, it is natural to identify S{f) as the scattering 
matrix in the case of arbitrary sources that give rise to the intensity of f. 

Based on QFT point of view, relation (PU]) indicates the appearance of the terms 
containing nonquantum fields that are characterized by the operators r{p^). Hence, 
there are terms with f in the matrix elements, and these f cannot be realized via 
real particles. The operator function r(p^) could be considered as the limit on an 
average value of some quantum operator (or even a set of operators) with an intensity 
that increases to infinity. The later statement can be visualized in the following 
mathematical representation: 



where a is the coherence (chaotic) function that gives the strength of the average 
S(p^,p^). 

In principal, interaction with the fields described by r is provided by the virtual 
particles, the propagation process of which is given by the potentials defined by the 
f operator function. 

The condition Mch (or Qo{R) ^ jh. * oo) in the representation 

ch 

lim S(p^,p')= lim 17o(i?) «((!>, /3) exp(-gV2) -^n(w,/3), 
with ^ 

means that the role of the arbitrary source characterized by the operator function 
f{pf_t) in b{pfj,) = a(p^) + f (p^) disappears. 



5 Green's function and kernel operator 

Let us go to the thermal field operator <I>(x) by means of the linear combination of 
the frequency parts t^"*" (x) and </>" {x) 

$(x) = -i= [0+(x)+0-(a;)] (24) 

composed of the operators 6(p^) and 6+(p^) as the solutions of equation and 
conjugate to it, respectively: 

^'(-)-/ (2.)32(p''L^)V^ ^"(^^)-"^ 
^"^(-)=/ (2.)32(p''Lt/- ^^"-^^""^- 
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The function ^(x) obeys the commutation relation 

[$(x),<i>(y)]_ =-zZ)(x) 

with [14] 



D{x)^:^e{x") 

Z TT 



Six^ 



e 



(a;2) Ji (my^^ 



where e(a;°) and 8(x^) are the standard unit and the step functions, respectively, 
while Ji{x) is the Bessel function. On the mass-shell, D{x) becomes [14] 



Dix) 



27r 



5ix')-^Qix') 



One can easily find two equations of motion for the Fourier transformed operators 
b{p^) and b+{p^jm S{^i) 

[uj - K{p^)]b{pf,) = F{p^) + p{ujp,e), (25) 

[lu - K+ip^)]b+ip^) - F+{p^) + p*{lop, e), (26) 
which are transformed into new equations for the frequency parts 0^(x) and (j)^{x) 



(27) 
(28) 



of the field operator ^{x) (|24]) 

ido^+ix) + I K{x- y) c^+(y)dy = f{x) 



idoc^-{x)+ K+{x-y)r{v)dy = f+{x), 



where 



and 



/(^) = 



d^p 



(27r)3(p2 +^2)1/ 
d^p 



— [F(p)+p{ujp,e)]t 



^"(^)-y (2.)3(,-.:^.)v. [^"w+^^(-^--)]-'^^- 

Here, the field components (t)'^{x) and 4)^{x) are under the effect of the nonlocal 
formfactors K{x — y) and K^{x — y), respectively. In general, these formfactors can 
admit the description of locality for nonlocal interactions. 

At this stage, it must be stressed that we have new generalized evolution Eqs. (I27p 
and (pS)) . which retain the general features of the propagating and interacting of the 
quantum fields with mass m that are in the heat bath (reservoir) and are chaotically 
distorted by other fields. For further analysis, let us rewrite the Eqs. P7)) and (|28p 
in the following form: 



ido(l)+{x) + K{x) * (j)+{x) = f{x), 
id^(t>- [x) +K+{x)*4>-{x)^ f+{x), 



(29) 
(30) 
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where A{x) ★ B{x) is the convoluted function of the generahzed functions A{x) and 
B{x). Applying the direct Fourier transformation to both sides of Eqs. (|29p and (|30p 
with the following properties of the Fourier transformation 

F[K{x) * = F[K{x)]F[^+{x)l 

we get two equations 

[p'> + Kip,M+ip,) = F[fix)l (31) 

[^p^ + K+{p,)]4>-{p,)^F[f+{x)]. (32) 

Multiplying Eqs. ([3T|) and (|43l) by — p*^ + A'+(pp) and + K{p^), respectively, we 
find 

i-p" + K+ip,)][p" + K{p,mP,) = Tip,), (33) 

where 

Tip,) = [V + K+ip,)]F[fix)] + [/ + Kip,)]F[f+ix)]. 
We are now at the stage of the main strategy: we have to identify the field $(x) 
introduced in Eq. (fT5|) and the field ^{x) built up of the fields 0+ and (j>~ as the 
solutions of generalized Eqs. (12711 and ([28]) . The next step is our requirement that 
Green's function G'(p^) in Eq. (|17p and the function r(p^), that satisfies Eq. 

hpO + i^+(pj]b° + i'(p^)]f(pj = 1, (34) 

must be equal to each other, where the full Green's function G{p^,g^,m^) 

G(]3^) ^ G(p ,g ,m ) ~ — (35) 

— p'^ — It 

has the same pole structure at p^ = m? as the free Green's function [14] with g being 
the scalar coupling constant and ^ is the one-loop correction of the scalar field. The 
dimensioneless function 1 — £,{p^ ^ 'm-'^) is finite at p^ = rn^. 

We define the operator kernel Kip,) in (P5)l from the condition of the nonlocal co- 
incidence of the Green's function G(p^) in Eq. (|17p . and the thermodynamic function 
Tip,) from jMl) in Si^^) 

F[Gip,) - tip,)] - 0. 
We can easily derive the kernel operator Kip,) in the form 

^ ^P> = ] 277^ 2^ (36) 

1 — g z,ip , m'^) 

where [14] 

e(m2) = 

and 

ap\^-) = -L_ 
The ultraviolet behaviour at ]p^ ] > > leads to 




f(p2,m2) 



4 



32 7r2 p2 



If 



,21 



ln^-^-.7re(p2) 
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6 Stochastic forces scale 



In paper [4] it has been emphasized that two different scale parameters are in the 
model which we consider here. One of them is the so-called "correlation radius" R 
introduced in (O and ® with © and HI]), ((H]). In fact, this i?-parameter 

gives the pure size of the particle emission source without the external distortion and 
interaction coming from other fields. The other (scale) parameter is the stochastic 
scale Lst which carries the dependence of the particle mass, the a-coherence degree 
and what is very important — the temperature T-dependence: 



1 



a{N)\p'^ ~ k{p)\ n{m,l3) 



(37) 



It turns out that this scale Lgt defines the range of stochastic forces acting the particles 
in the emission source. This effect is given by a(iV)-coherence degree which can be 
estimated from the experiment within the two-particle BE correlation function C2 (Q) 
as Q close to zero, C2(0), at fixed value of mean multiplicity (N): 



a{N) 



2-^2(0) + V2-C2(0) 



, C2(0) = C2(0)/e(iV). 



C2(0)-l 

In formula p7p n{m, (3) is the thermal relativistic particle number density 



n(?n, /3) = 3 



(2^ 



n{uj, P) — 3- 



2 7r2 



1=1 



(38) 



(39) 



where K2 is the modified Bessel function. For definite calculations we consider cor- 
relations between charge pions. The result can be extended to heavy particles case, 
e.g., for charge and neutral gauge bosons that is essential program for the LHC. 
The stochastic scale Lst tends to infinity in case of particles are on mass-shell, i.e., 
\p^~K{p)\ ~~f which enters the L^jS denominator (P7|). However, Lst will be bounded 
due to stochastic forces acting the particles where 



p" 



as 5^^(^2,771^) < 1. 

Within our aim to explore the correlation between charged pions, Lst has the form 



3a(7V)Ae2(^2+„,2)3/4^^^3/2 I ^ 



(40) 



where the condition I j3^J w? + fj? > 1 for any integer I in (|39p was taken into account. 
The only lower temperatures will drive Lst within the formula ()40p even if /i = and 
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I = 1 with the condition T < m. Note that the condition /i < m is a general 
restriction in the relativistic "Bose-Hke gas", and fi ~ m corresponds to the Bose- 
Einstein condensation. 

For large enough T no the dependence ol the chemical potential /i is found for Lgt ■ 



3C(3)a(7V) Ae2r3 



(41) 



where the condition T > /y/i^ + m^, I — 1, 2, ... is taken into account. The origin of 
formula ([4T|) comes from 



n{m, (3) — > n(/3) 



^C(3) 



(42) 



where neither a pion mass m- nor /i- dependence occurred; ^(3) = X^j^i — 1.202 
is the zeta- function with the argument 3. For high momentum pions {p^ ~ '^m?) the 
actual mass-dependence occurred for Lgt- 



L, 



at low T, and 



3a(N) m2 (^2 + ^2)3/4 ^t.) 



3/2 
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3C(3)a(iV) m2 



(43) 



(44) 



at high temperatures, if ^^^(to^) << 1, ^(to^) ^ ©(O.Ol/m^) and {fp' /Am?) « 1 are 



valid in both temperature regime cases. Formula (|42p reproduces the ~ behavior 
which is the same as the thermal distribution (in terms of density) for a gas of free 
relativistic massless particles. Such a behavior is expected anyway in high temperature 
limit if the particles can be considered as asymptotically free in that regime. 

Actually, the increasing of T leads to squeezing of Lgt, and Lst{T = Tq) = i? at 
some effective temperature Tq. The higher temperatures, T > Tq, satisfy to more 
squeezing effect and at the critical temperature Tc the scale Lst{T = Tc) takes its 
minimal value. Obviously Tc defines the phase transition where the deconfinement 
will occur. Since all the masses tend to zero (chiral symmetry restoration) and a ^ 
at T > Tc one should expect the sharp expansion of the region with Lst(T > Tc) ~^ oo. 
The following condition him, 13)-Vt^ = 1 provides the phase transition (transition from 
hadronizing phase to deconfinement one) with the volume ~ (4 7rr3/3), where 
is the pion charge radius. Actually, the temperature of phase transition essentially 
depends on the charge (vector) radius of the pion which is a fundamental quantity in 
hadron physics. A recent review on r^r values is presented in [15]. 

What we know about the source size estimation from experiments? DELPHI 
and L3 collaborations at LEP established that the correlation radius R decreases 
with transverse pion mass mt as i? ~ a + b/^/ml for all directions in the Longi- 
tudinal Center of Mass System (LCMS). ZEUS collaboration at HERA did not ob- 
serve the essential difference between the values of i? - parameter in tt^tt*, K^K^ 
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and K^K"^ pairs, namely = 0.666 ± 0.009{stat) + 0.022 - 0.033{syst.)fm, 

Rk,k, = 0.61 ± Om{stat) + 0.07 - 0.08(sysi.)/m and Rk±k± = 0.57 ± 0.m{stat) + 
0.15 — Q.06{syst.)fm, respectively. The ZEUS data are in good agreement with the 
LEP for radius R. However, no evidence for y/s dependence of R is found. It is 
evidently that more experimental data are appreciated. However, the comparison 
between experiments is difficult mainly due to reference samples used and the Monte 
Carlo corrections. 

Finally, our theoretical results first predict the Lgt in (|40p and (|¥T|) . and both 
mass- and temperature - dependence are obtained clearly. This can serve as a good 
approximation to explain the LEP, Tevatron and ZEUS (HERA) experimental data. 
We need that the pion energies at the colliders are sufficient to carry these studies out 
(since the Acp dependence). Careful simulation of their (pions) signal and background 
are needed. The more precisely measured pion momentum may be of some help. Also, 
determination of the final state interactions may clarify what is happening. 

7 Conclusions 

To summarize: we find the time dependence of the correlation function C2{Q) calcu- 
lated in time-dependent external field provided by the operator r(p, t) and the chaotic 
coherence function a{m,/3). Based on this approach we emphasize the explanation 
of the dynamic origin of the coherence in BEC, the origin of the specific shape of 
the correlation C2(Q) functions, and finding the dependence on the particle energy 
(and the mass) due to coherence function a, as seen from the QFTp. Actually, the 
stochastic scale Lst decreases with the particle energy (the mass m). It is already 
confirmed by the data of LEP, Tevatron and HERA (ZEUS) with respect to the size 
of particle source. 

In the framework of QFTg^ the numerical analysis of experimental data can be 
carried out with a result where important parameters of C2 (Q) and C2~^^{Q) 
functions are retrieved (e.g., Co, i?. A, e, A'^, a, Lst, T). 

The correlations of non-identical particles pairs can be observed and the corre- 
sponding C2 parameters is retrieved. The off-correlation effect is given by the space- 
time distribution (|14p containing the sum pi + P2, and this effect is sufficient if the 
factor containing the sum pi + p^ in (fTO|) is not too small. The off-correlation effect 
is possible if the particle energies p^{i = 1,2) are small enough. 

Besides the fact that like, e.g., tt^tt^ BEC the correlations tt^tt^ can serve as 
tools in the determination of parameters of the particle source. And besides the fact 
that these correlations play a particularly important role in the detection of random 
chaotic correction to BEC. 

The stochastic scale Lst decreases with increasing temperatures slowly at low 
temperatures, and it decreases rather abruptly when the critical temperature is ap- 
proached. 

We claim that the experimental measuring of R (in fm) can provide the precise 
estimation of the effective temperature Tg which is the main thermal character in the 
particle's pair emitter source (given by the effective dimension R) with the particle 

mass and its energy at given a fixed by C2{Q = 0) and {N). Actually, Tq is the 
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true temperature in the region of multiparticle production with dimension R = L^t , 
because at this temperature it is exactly the creation of two particles occurred, and 
these particles obey the criterion of BEC. 

We have found the squeezing of the particle source due to decreasing of the corre- 
lation radius R in the case of opposite charge particles. The off-correlated system of 
non-identical particles is less sensitive to the random force influence (a- dependence) . 

The results obtained in this paper can be compared with the static correlation 
function (see, e.g., [16] and the references therein relevant to heavy ion collisions). 

Finally, we should stress a new features of particle-antiparticle BEC which can 
emerge from the data. It is a highly rewarding task to experimental measurement of 
non-identical particles. 

There is much to be done for C2{Q) investigation at hadron colliders. The time 
is ripe for dedicated searches for new effects in C2{Q) function at hadron colliders to 
discover, or rule out, in particular, the a{N) dependence. 

In conclusion, the correlations of two bosons in 4-momentum space presented in 
this paper offer useful and instructive complimentary viewpoints to theoretical and 
experimental works in multiparticle femtoscopy and interferometry measurements at 
hadron colliders. 
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